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ABSTRACT 


Temporal  voltage-current  dynamic  characteristics  have  been  com¬ 
puted  using  the  Townsend  continuity  equations  and  Poisson's  equation, 
with  suitable  initial  and  boundary  conditions.  In  the  current  region 
of  the  negative  static  characteristic,  small  current  perturbations 
initially  grow  exponentially  in  time  if  the  external-circuit  time  con¬ 
stant  is  greater  than  the  tube  current -growth  time  constant.  In  a 
limited  range  of  applied  voltages,  these  oscillations  become  precisely 
repetitive  in  time,  corresponding  to  the  well  known  relaxation  oscil¬ 
lations.  The  frequencies  of  the  oscillations  obtained  from  these 
calculations  agree  with  those  measured  at  HDL  more  closely  than  do 
those  obtained  from  the  traditional  published  formula.  The  inductances 
and  the  negative  resistances  calculated  from  small -amplitude  decaying 
oscillations  agree  precisely  with  those  calculated  using  a  sinusoidal 
incremental  applied  voltage  of  the  same  frequency  and  average  current. 
Qualitative  agreement  is  also  obtained  with  laboratory  measurements  of 
these  quantities.  Relaxation  oscillations  are  also  calculated  psing 
voltages  that  closely  approximate  the  peaks  of  60-cps  applied  voltages. 
Most  of  the  calculations  presented  are  for  an  argon  gap  of  1  cm  and 
for  a  pressure  of  2  torr.  Experimental  data  are  presented  for  both 
argon  and  hydrogen. 

1.  INTRODUCTION 

The  well  known  phenomenon  of  relaxation  oscillations  in  cold- 
cathode  gas  tubes  is  of  great  practical  value  (ref  1) .  In  most  treat¬ 
ments  of  these  oscillations,  the  gas  tube  is  considered  as  a  perfect 
switch,  which  closes  to  discharge  a  condenser  from  the  tube  breakdown 
voltage  to  its  extinction  voltage  and  then  opens  in  zero  time.  This 
simplified  theory  may  lead  to  large  errors  in  calculating  the  frequency 
of  the  relaxation  oscillations  (ref  2).  This  report  gives  a  fuller 
treatment  of  a  paper  presented  at  the  1962  summer  meeting  of  the 
American  Physical  Society  (ref  3),  which  for  the  first  time  presented 
a  rigorous  calculation  of  the  processes  within  the  gas  tube  in  a  relaxa¬ 
tion  oscillator  circuit. 

2.  CALCULATIONS 

2.1  Model 


The  model  and  resulting  formulation  used  in  the  calculations, 
now  adapted  for  the  IBM  7090,  have  been  presented  previously  (ref  4,5). 

In  brief,  the  Townsend  continuity  equations  for  gas  ionization  are 
solved  simultaneously  with  Poisson's  equation  to  account  for  space  charge 
effects.  The  computer  output  includes  the  variation  of  the  tube  current 
and  voltage  with  time.  The  boundary  conditions  on  the  gas  tube  are  im¬ 
posed  by  a  series  resistance  R  and  a  parallel  capacitor  C,  with  a  volt-  . 
age  supply  which  permits  either  a  step  voltage  or  an  incremental  sinusoidal 
voltage  to  be  added  to  the  constant  d-c  supply.  The  initial  conditions 
are  arbitrary,  but  a  stationary  state  is  generally  used. 
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2.2  Approximations 


The  frequency  of  relaxation  oscillations  has  been  routinely 
calculated  as  the  reciprocal  of  the  time  required  for  the  condenser 
to  charge  through  the  series  resistance  from  the  extinction  voltage 
of  the  tube  to  the  breakdown  voltage  Vb>  This  time  is  readily 
calculated  as 

T  «  h  C  In  [(V  -  V  )/(V  -  V.  )]  (1) 

s  a  x  a  b  4 

where  V  is  the  applied  voltage.  This  equation  gives  no  inf onset ion 
on  the  range  of  values  of  Vg,  Rg,  or  C  for  which  oscillations  occur, 

or  on  the  currents  involved. 

A  criterion  for  the  stability  of  the  gas  tube  has  been 
given  by  Kaufmann  (ref  6).  The  following  analysis  was  made  independ¬ 
ently  but  is  similar  to  that  of  Kaufmann.  It  has  been  found  (ref  5) 
that  the  impedance  properties  of  a  gas  discharge  may  be  approximated 
by  an  inductance,  L,  and  a  differential  resistance,  R,  The  effective 
circuit  is  shown  in  figure  1. 


dVT 

diT 


Figure  1.  Effective  circuit  used  to  analyze  gas 
tube  oscillations i  The  tube  is  rep¬ 
resented  as  an  Inductance  and  an  a-c 
resistance. 
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Let  the  steady-state  tube  current  be  I  and  the  current  at  any  time 

be  I(t) .  Then  the  differential  equation  for  the  difference  current, 
i(t)  a  I(t)  -  is  given  by 


(2) 


The  derivation  of  this  equation  and  its  solution  is  given  in  the 
appendix.  One  form  of  the  solution  is 

i(t)  a  *(<?)exp(Xt)  sin(u)t  +$)  (3) 

where  X  =  -  \  g  +  ~-]  and  &  -  (l  +  -f-}  -  X» 

s  s' 

From  equations  (2)  and  (3)  it  is  seen  that  if  R,  the  slope  of  the 
static  characteristic,  is  positive,  the  tube  is  stable  for  small 
perturbations.  However,  if  R  is  negative  and  of  greater  magnitude 
than  R  ,  the  tube  is  unstable.  This  Instability  corresponds  to  the 
breakdown  of  the  tube  into  the  glow  discharge  mode.  If  R  is  negative 
and:  (1)  if  the  internal  time  constant  of  the  tube  L/R  equals  in 
magnitude  the  external  circuit  time  constant  R^C,  then  X  *  0  and  the 

tube  oscillates  with  a  constant  amplitude,  with  the  frequency  given 
by  the  angular  frequency,  U) ;  (2)  if  L/-R  is  greater  than  RgC,  then 

any  perturbation  in  the  tube  current  dies  out  exponentially;  (3)  if 
L/rR:iJ!  less  than  R  C,  then  X,  the  current  growth  constant,  is  positive 
and  any  perturbation  will  initially  grow  exponentially  in  time.  How¬ 
ever,  both  R  and  L  are  functions  of  the  tube  current,  R  becoming  large 
and  positive  in  the  abnormal  glow. 

Equations  (3)  may  be  solved  for  R  and  L  to  obtain 


CR 


L  =  - 


1  +  2  R  CX 
s 


R  = 


-R  (1  +  2  R  GX) 
s  s 


1+2  RSCX  +  I^C®  <00®  +  xa) 


(4) 


3.  HrraiMKNTAL  MEASUREMENTS 

Experimental  measurements  for  relaxation  oscillations  have  been 
reported  in  reference  2.  The  present  paper  reports  more  extensive 
work  and  analysis  completed  at  HDL  since  that  report.  The  experimental 
arrangement  is  shown  in  figure  2  (next  page) •  Data  were  obtained  using 
two  different  experimental  parallel-plkte  tubes.  One  tube  contained 
argon  gas  at  a  pressure  of  9  torr  with  nickel  electrodes  spaced  0.84 
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cm  apart.  The  second  tube  was  built  with  a  nickel  cathode,  with  the 
possibility  of  varying  both  the  filling  gas  pressure  and .the  electrode 
spacing. 
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Figure  2.  Schematic  circuit  of  experimental 
apparatus  used  to  measure  tube 
oscillations. 


Relaxation- type  oscillations  were  observed  and  studied  in  the 
course  of  obtaining  the  static  characteristic  curve  in  a  current  range 
from  10  *to  10"3  amp,  which  includes  the. negative  slope  region  of 
the  curve . 

I  4.  RESULTS 


4.1  Calculations 


Most  of  the  calculations  have  been  made  for  argon  at  a 
pressure  p  of  2  torr,  a  gap  distance  d  of  1  cm,  and  a  secondary  ioniza¬ 
tion  coefficient  y  equal  to  0.04.  The  variation  of  the  primary  ioniza¬ 
tion  coefficient,  the  ion  and  electron  mobility  with  the  field  E^ is 
given  in  reference  4. 


) 


An  example  of  the  buildup  of  oscillations  from  a  small 
perturbation  is  shown  in  figure  3.  After  a  varying  number  of  oscil¬ 
lations,  the  oscillations  become  precisely  repetitive  in  time. 

Using  Rg  =  10®  ohms,  C  =  3  x  10“®  f,  a  discharge  area  is  =  10  cm®,  and 
an  externally  produced  current,  J  =  2  x  10~®amp/cm,  relaxation  oscil¬ 
lations  were  found  to  occur  witli  applied  voltages  of  173.5  to  180.5  v. 

At  both  173.4  and  181  v  the  oscillations  die  out..  To  define  the 
boundary  values  more  closely  would  take  excessive  computer  time.  The 
variation  of  both  tube  current  and  voltage  with  time  is  shown  in  figure 
4,  while  a  plot  of  tube  voltage  versus  current  is  given  in  figure  5. 

It  will  be  noted  that  the  frequency  of  oscillations  increases  with  ap¬ 
plied  voltage.  The  variation  of  frequency  with  applied  voltage  is 
given  in  figure  6  (a),  together  with  that  predicated  by  equation  (1). 

It  is  seen  that  the  frequency  calculated  from  equation  (1)  is  about 
5  to  6  times  that  calculated  from  the^cbmputet  data.  Measurements  re¬ 
ported  in  reference  2  indicated  that  equation  (1)  gives  frequencies 
about  3  times  greater  than  those  measured  under  somewhat  different 
experimental  conditions.  A  plot  of  frequency  versus  the  time  average 
of  tube  current  is  given  in  figure  6  (b),  on  a  log-log  plot.  The 
points  corresponding  to  the  relaxation  oscillations  drop  slightly 
below  the  straight  line  formed  by  points  obtained  outside  this  region. 

If  JQ  is  reduced,  the  frequency  of  the  lower  applied  voltage  oscilla¬ 
tions  drops  sharply,  requiring  excessive  computer  time.  At  higher 
applied  voltages,  the  current  does  not  reach  as  low  a  minimum  value 
and  the  J  value  has  less  effect.  In  fact  this  higher  minimum  current 
is  effective  in  lowering  the  breakdown  voltage  in  the  same  way  as  is 
J  .  Since  the  tube  breaks  down  at  a  lower  potential,  less  charge  is 
scored  in  the  capacitor  and  consequently  the  current  peak  is  not  so 
large.  This  may  be  noted  in  figures  3  and  4.  The  static  character¬ 
istic  calculated  from  runs  of  other  types  is  shown  in  figure  5.  It  is 
seen  that  the  static  characteristic  separates  the  region  where  current 
increases  with  time  and  that  where  current  decreases  with  time.  It  may 
also  be  noted  that  the  load  line  crosses  the  dynamic  paths  at  the  maximum 
and  the  minimum  voltages.  The  time  average  currents  and  voltages  of  the 
relaxation  oscillations  are  plotted  as  individual  points  in  figure  5. 

The  buildup  or  decay  times  of  the  relaxation- type  oscillations 
and  their  frequency  may  be  used  to  compute  the  tube  Inductance  L,  and 
a-c  resistance  R,  by  means  of  equation  (4) .  An  alternative  method  of 
determining  L  and  R  by  forced  oscillations  of  an  applied  frequency  F, 
was  presented  in  reference  4.  The  good  agreement  of  the  Inductances 
calculated  by  the  two  methods  is  presented  in  figure  7,  where  the  LI 
product  is  plotced  against  the  current  density.  The  frequency  for  the 
forced  oscillations  is  1  x  104  cps,  while  the  self-oscillations  vary 
somewhat  from  that  figure.  The  "no  space  charge  approximation"  is  dis¬ 
cussed  in  reference  4.  The  minimal  value  of  the  LI  product  is  attained 
at  a  current  density  about  125  percent  of  the  normal  current  density. 
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Figure  3.  Buildup  o*f  oscillations  from  a  small  perturbation  calculated  for  an  argon  tube 
operating  in  the  negative  resistance  region. 
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(a)  Calculated  variation  of  frequency  with  applied  voltage. 
Curve  marked  "RC"  la  obtained  from  equation  1.  Results 
obtained  for  stable  relaxation  oscillations  are  shown  by 
open  point b,  while  decaying  oscillation  frequencies  are 
shown  by  crosses,  ,  . 
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(b)  Calculated  variation  of  frequency  against  average  current. 
Open  points  represent  stable  oscillations  while  crosses 
Indicate  decaying  oscillations. 


Figure  6.  Frequency  as  a  function  of  applied  voltage  and  current. 


The  variation  of  impedance  and  the  phase  angle  (phase  lag  of 
the  current  from  the  voltage)  with  frequency  is  shown  in  figure  8  for 
two  current  densities,  14 .0  ua  cm*  (approximately  the  normal  current 
density)  and  17.4  |ia/cm8 .  It  is  noted  that  phase  angle  reaches  a 
maximum  at  roughly  half  the  frequency  at  which  the  impedance  attains  its 
maximum. 


The  variation  of  the  LI  product  with  frequency  is  shown  in 
figure  9  for  currents  in  the  subnormal  glow,  the  normal  glow,  and  the 
abnormal  glow.  Note  that  the  maximum  shifts  to  higher  frequencies  as 
tube  current  increases.  A  plot  of  U)LIp  versus  Rip  is  given  in  figure  10. 
Experimental  plots  of  d)L  versus  R  have  generally  given  semicircular 
curves.  Figure  10  shows  that  the  similarity  laws  are  obeyed  over  the 
limited  range  of  current  and  pressure  tested . 

It  has  been  found  that  the  portion  of  the  60-cps  a-c  voltage 
above  the  breakdown  voltage  may  be  closely  approximated  by  calculations 
using  an  approach  voltage  and  an  Incremental  sinusoidal  voltage  of  ap¬ 
propriate  amplitude  and  frequency.  Figures  11  and  12  show  the  results 
of  a  sequence  of  such  calculations  simulating  60-cps  a-c  voltages  of 
the  various  amplitudes  noted.  Aside  from  the  applied  voltage,  the  cir¬ 
cuit  and  tube  parameters  are  the  same  as  those  used  in  the  previous 
figures.  The  top  portions  of  each  figure  show  the  applied  voltage  and 
the  tube  voltage  as  a  function  of  time.  The  middle  figures  present  the 
tube  and  external  resistor  current  as  a  function  of  time.  The  lower 
figures  give  the  dynamic  voltage-current  path  and  the  static  character¬ 
istic.  It  should  be  noted  that  the  number  of  pulses  increases  with 
increasing  voltage,  that  a  pulse  is  nearly  full  sise  if  it  exceeds  the 
"breakdown  current,"  and  that  successive  current  pulses  are  smaller 
than  the  initial  breakdown  pulse.  Closer  scrutiny  reveals  that  the 
frequency  is  greater  at  higher  applied  voltages.  When  the  load  line 
crosses  the  static  characteristic  in  the  abnormal  glow  region,  the 
oscillations  rapidly  die  out.  The  current  dies  out  without  oscilla¬ 
tions  building  up  as  the  applied  voltage  decreases  to  below  the  break¬ 
down  value.  At  the  two  highest  applied  voltages,  calculations  were 
extended  only  to  approximately  the  maximum  voltages,  since  otherwise 
the  computer  time  would  be  excessive. 

4 .2  Experimental 


Current-voltage  characteristics  were  obtained  using  d-c  meters. 
Therefore,  when  oscillations  are  observed  on  the  oscilloscope,  the  cur~ 
rent-voltage  points  must  be  considered  as  representing  "average"  read¬ 
ings.  Such  a  characteristic  is  shown  in  figure  13,  where  the  oscillatory 
regions  are  noted  on  the  graph. 

Figure  14  shows  the  variation  of  the  frequency  of  relaxation 
oscillations  in  argon  as  a  function  of  the  "average"  current.  The 
data  were  taken  at  a  pressure  of  9  torr  and  a  spacing  of  0.84  cm.  The 
effect  of  the  circuit  shunt  capacitance  on  the  frequency  of  oscillation 
is  illustrated  here.  One  of  these  curves  has  been  replotted  in  figure 
15  against  applied  voltage,  and  compared  with  that  curve  predicted  by 
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Figure  10%  Calculated  variation  of  reactahce-current^pxeeeure  product  with 
reaiatance-current-preaaure  product.  Raraaeter  4»  the  applied 
frequency. 
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Figure  14,  Measured  variation  of  frequency  of  relaxation  oscillations 
as  a  function  of  the  average  current  for  a  series  resist¬ 
ance  of  10  meghoms  and  various  values  of  shunt  capacitance 
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Figure  15.  Variation  of  frequency  of  oscillations  with  applied  voltage. 

Curve  labeled  "RCrv  was  calculated  using  equation  1,  while 
•  points  show  measured  values. 


equation  (1).  Note  that  the  frequency  of  the  predicted  "R C"  curve  is 
3  to  4  tines  that  of  the  actual  experimental  data. 

Similar  measurements  made  in  hydrogen  with  a  pressure  of  2,2 
torr  and  a  gap  distance  of  1  cm  are  shown  in  figure  16.  This  graph 
plots  frequency  against  current  and  shows  frequency  as  a  function  of 
both  series  resistance  and  shunting  capacitance. 

To  correlate  the  data  taken  in  argon  with  that  in  hydrogen, 
the  information  from  the  previous  figures  has  been  replotted  together 
in  figure  17,  where  the  frequency-pressure  product  has  been  plotted 
against  current . 

Experimental  work  has  also  been  done  with  a  60-cps  voltage 
applied  to  the  tube  to  compare  with  similar  computer  calculations 
reported  above.  Examples  of  oscilloscope  traces  of  this  work  are 
photographed  in  figure  18.  This  series  of  pictures  was  taken  with 
argon  at  a  pressure  of  9  torr  and  a  gap  distance  of  0.84  cm. 

Experimental  work  on  the  impedance  porperties  of  gas  dis¬ 
charges  has  been  reported  previously  (ref  5) . 

5.  DISCUSSION 

The  traditional  model  for  calculating  the  frequency  of  relaxation 
oscillations  in  gas  tube  circuits  assumes  that  the  tube  breaks  down  at 
a  fixed  breakdown  voltage  and  extinguishes  at  a  fixed  voltage  (the 
normal  glow  voltage)  in  a  time  negligible  compared  with  the  time  for 
the  condenser  to  recharge.  The  frequency  may  then  be  calculated  by 
equation  (1).  It  was  shown  in  reference  2,  and  is  further  substan¬ 
tiated  in  this  report,  that  the  frequencies  thus  calculated  are  larger, 
by  a  factor  of  about  3  or  4,  than  those  measured  experimentally  at  low 
pressures.  An  examination  of  figures  4  and  5  shows  that  the  error  may 
stem  from  three  sources.  Pirst,  the  breakdown  voltage  is  not  constant 
but  depends  upon  both  the  approach  current  and  the  rate  of  voltage 
rise.  It  is  normally  considerably  higher  than  the  breakdown  voltage 
as  measured  quasi-statically .  Likewise,  the  extinction  voltage  is  not 
constant  but  is  considerably  below  the  normal  glow  voltage.  Normally, 
the  higher  the  breakdown  voltage,  the  larger  the  current  spike  and  the 
lower  the  extinction  voltage.  Finally,  the  time  for  the  breakdown 
cannot  be  neglected  at  low  pressure. 

Certain  limitations  of  the  calculations  and  of  the  ability  to 
make  measurements  under  ideal  conditions  prohibit  direct  comparison 
of  the  two  results.  First,  the  calculations  are  one-dimensional 
while  the  discharge  is  not.  Some  success  has  been  achieved  in 
measuring  cross-sectional  area  in  static  discharges,  and  this  informa¬ 
tion  is  of  some  value  in  interpreting  dynamic  measurements.  Second, 
the  calculations  are,  at  present,  limited  to  a  constant  secondary 
ionization  coefficient  y ,  Independent  of  the  cathode  field  and  of  the 
current  flowing.  It  is  generally  accepted  that  these  factors  cause 
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I  -  AMPERES 


Figure  16,  Measured  variation  of  frequency  with  average  current  in 

Hggap  at  a  pressure  of  2.2  torr.  Series  resistance  values 
in  megohms  are  given  by  the  following  symbols:  crosses  -1 
circles  -  10,  triangles  4,  and  plus  signs  -  50. 


4  6  a  10“ 8  2  4  e  8|0“*  2 
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Figure  17.  Measured  frequency-pressure  product  as  a  function  of 
current  for  both  argon  and  hydrogen. 


Figure  18  .  Photographs  of  relaxation  oscillations  occurring  with  60  cps 
applied  voltage.  Both  current  and  voltage  traces  are  shown 
in  each  photograph.  Successively  increased  number  of  oscil¬ 
lations  are  observed  at  increasing  applied  voltages. 


large  changes  in  y  at  low  pressures.  Moreover,  the  calculations  have 
assumed  that  v  is  entirely  due  to  the  effect  of  ions  striking  the  cathode, 
with  no  allowance  for  a  photon  effect.  Since  this  is  the  case,  all  cal¬ 
culations  reported  herein  have  completely  ignored  the  transit  time  of  the 
electrons,  for  reasons  outlined  elsewhere  (ref  8). 

Despite  the  above  limitations,  the  agreement  between  the  experimental 
data  and  the  calculations  is  quite  promising.  As  noted  above,  the  cal¬ 
culated  frequencies  of  relaxation  oscillations  are  about  one-sixth  of 
that  calculated  using  equation  (1)  with  the  calculation  parameters,  while 
the  experimental  measured  values  are  about  one-third  of  ttypse  calculated 
using  equation  (1)  under  the  experimental  conditions.  This  indicates 
an  error  of  about  a  factor  of  about  two  between  the  calculations  and  the 
experimental  data.  A  few  exploratory  calculations  have  shown  that  an  as¬ 
sumption  of  50  percent  contributions  to  the  secondary  ionization  coef¬ 
ficient  from  each  of  the  ion  component  and  the  photon  component  is  able 
to  cut  this  error  in  about  half.  Further  calculations  will  be  necessary 
to  verify  this. 

A  more  direct  comparison  can  be  made  between  the  calculated  and 
the  experimental  variation  of  frequency  versus  current.  To  make  this 
comparison  it  is  profitable  to  use  similarity  parameters,  frequency/pres¬ 
sure,  F/p,  and  current/( pressure)3 ,  I/p? .  Such  a  plot  is  shown  in  figure 
19  for  the  calculated  data  shown  in  figure  6  (b)  and  the  experimental 
data  for  one  R  C  value  taken  from  figure  17.  The  agreement  of  both  the 
slope  and  intercept  is  satisfactory. 

The  variation  of  the  LI  product  as  a  function  of  frequency,  as 
shown  in  figure  9,  shows  surprising  agreement  with  the  measurements  of 
Crawford  (ref  9)  in  hot  cathode  mercury  discharges  at  1-p,  pressure.  ' 
Figure  3  of  reference  9  shows  that  the  measured  LI  product  is  1;2  x  10"4 
from  10  to  30  kc,  then  drops  nearly  linearly  on  a  log-log  plot  to  1  x  10~6 
at  1.8  Me.  This  agreement  is  remarkable  considering  the  dissimilarity 
of  the  gas,  pressure,  cathode  mechanism,  etc.  The  semicircular  impedance 
locus,  similar  to  that  calculated  and  shown  in  figure  10,  was  also  noted 
by  Crawford  (ref  9). 

The  voltage  and  current  variation  with  time, calculated  to  simulate, 
the  60-cps  discharge  and  shown  in  figures  11  and  12,  shows  good  agreement 
with  the  experimental  data  shown  in  figure  18  under  comparable  discharge 
conditions.  The  experimental  data  of  El-Bakkal  and  Loeb  (ref  10)  taken 
in  argon  in  glass  cells  also  show  complete  qualitative  agreement  with 
the  calculations  of  figures  11  and  12. 

6.  CONCLUSIONS  AND  FUTURE  PROGRAM 

It  has  been  shown  that  the  dynamic  Townsend  equations  for  continuity 
in  a  gas  discharge,  with  suitable  boundary  conditions  and  including 
space  charge  effects  are  able  to  explain  relaxation  oscillations  in  gas 
tube  circuits.  Although  direct  comparison  with  experimental  measurements 
has  not  yet  been  accomplished,  qualitative  agreement  is  excellent  and 
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Figure  19.  Comparison  of  computed  frequency/ pressure  quotient  variation  with 
current-pressure  squared  quotient  with  measured  values.  Com¬ 
puter  results  for  argon  are  shown  as  triangles,  measured  results 
in  argon  by  crosses, and  in  hydrogen  by  circles. 


the  variation  of  frequency  with  current  la  in  quantitative  agreement . 

It  is  planned  to  make  further  calculations  to  determine  the  effect 
of  the  variation  of  pressure,  gap  length,  ionisation  coefficients,  and 
a  number  of  parameters  peculiar  to  the  computer  program,  upon  the  cal¬ 
culated  relaxation  oscillations.  Further  experimental  work  will  be 
done  whenever  it  seems  feasible  to  check  the  calculations. 
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AHMIXX— Dill  VAT  I  ON  AND  SOLUTIONSOT  DIFFERENCE  CURRENT  EQUATION 

The  fallowing  definitions  are  made  (see  figure  1)  : 

C  :  stray  capacitance  plus  plate  capacitance, 
effective  Inductance  of  the  tube. 

Rt :  instantaneous  tube  resistance  (slope  of  the  static  characteristic). 

R  :  series  resistance, 
a 

V 

to:  maintaining  voltage. 

V^:  voltage  across  the  tube. 

*T:  current  through  the  tube. 

*1:  current  through  Rs* 

Ig :  current  through  C. 

:  stationary  current  through  tube. 
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Differentiating  equation  (5)  with  respect  to  time  t,  we  see  that 


•i-'s 


Substituting  (7)  into  (6) 

d2i_  di 

W  ^2“  *  +  W  —  +  <Rs  +  V1!  =  Vo  -  <Rs  +  V1., 


d  iT  .  S  +  Vt  diT  Rs  +Rt 


dt  r8ltc 


dt  R  L  C 
ST 


v  ■  v«  -  *  V1. 


=  0  by  definition  of  I  . 


This  is  the  differential  equation  that  must  be  solved.  The  fact  that 
the  equation  has  constant  <poeff icients  makes  the  solution  and  analysis 
much  easier. 

First,  for  convenience,  let  us  make  the  following  definitions  and 
substitutions: 

X  -  *  j-i.)  .nd  .2  =  (1  •  ♦  jig  -  X2  (9) 

T  s  s  ^T 

Now  the  differential  equation  (8)  becomes 


d2 i_  di_ 

- 1  -  2X  — ^  +  (u)2  +  X2)i  =  0 

*+*  dt  T 


The  general  solution  to  (10)  is 

iT  =  A  exp{X-(-U)2)1/2}t  +  B  exp{X  +  (-tt)2)1^2^ 


where  A  and  B  are  arbitrary  constants.  To  determine  the  character¬ 
istics  of  expression  (11),  it  is  necessary  to  examine  the  solution 
rather  closely,  keeping  in  mind  the  many  possible  values  of  the  circuit 
components  of  figure  1. 
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Case  L  R>  0 


From  the  definition  of  X  in  expression  (9),  we  see  that  if  R^ 
a  positive  number ,  then  X  is  always  negative, 

A.  u>2  negative  {X2  > 


Under  this  condition,  the  exponentials  in  solution  (11)  are  real. 
Therefore^  the  solution  takes  the  form 

(i_  =  A  exp  at  n  =  1,2), 

T  n  n  7 

2  1  /2  2  1  /2 

where  a^  =  X  -  )  and  ag  =  X  +  (-0)  )  ,  a^^  and  ag  always  being 

negative  numbers. 


B.  0)  a  0 


{x2  =  (1  +  ?)  "} 


R  '  L  CJ 
s  T 


In  this  case  there  is  a  double  root  of  the  characteristic  equa¬ 
tion,  that  root  being  X,  Therefore,  the  general  solution  becomes 


(iT  =  (A  +  Bt)  exp  Xt) 


C,  a)2  positive  [X2  <  (1  +  — )  1—) 


’  r  ^  R  L,C 

S  T 

Here  the  exponentials  in  (11)  are  imaginary.  One  form  of  the 
general  solution  is  then 

(iT  a  A  exp  (Xt)  sin  (u>t  +  B) 

This  is  an  oscillating  function  in  time  with  expontially  decreasing 
amplitude. 

Case  II*  R^=s0 

In  this  case,  X  again  is  always  negative,  and  the  forms  of  the 
solution  are  the  same  as  in  Case  I,  except  that  the  definitions  of  X 
and  u)  become 


.  1  ^  2  1  2 

X  -  -  2RC  and  »  "  IJ-  ~  X 

S  T 
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Case  III.  0 

A.  X  =  0  (ly/L, 


T'  R  C 
s 


The  differential  equation  (10)  becomes 


d2i 


dt 


T  2 

4  +  U)i  =  0 
2  T 


2 

1*  OJ  negative  {  | |  >Rg} 

The  solution  is  now  real,  taking  the  form 


(iT  =  An  ®*P  ^ 


2.  u)2  positive  {[R^l  <  Rg} 


The  exponentials  in  (11)  are  imaginary,  one  form  of  the 


solution  being 

(iT  =  A  sin  (u)t  +  B) 

a  constant-amplitude  oscillation  with  frequency 


R  +RA1/2 


/VRt\ 

Iwl 


B.  X  <  0 


1  1  Lj  RCJ 

T  S 


This  case  has  the  same  form  as  Case  I,  since  X  <  0, 
C.  X  >  0 


fl-cl  >Fc} 

T  s 


1.  U)2  =  0 


-  (i  ♦  re ■< 


Again  the  characteristic  equation  has  a  double  root  (X),  and  the  solution 
is 


(i  =  (A  +  Bt)  exp  Xt> 


2.  u)  positive  \ 


1*4.1 


<  R 


7  and 


V 


>  X 


35 


=  A  exp  Xt  #in  («t  + 

>■  -  i r 

an  oscillating  current,  with  an  exponentially  rising  amplitude. 
3.  u)2  negative 


U)2  is  negative  (1)  if  | |  >  Rg  orf2)  if  (r^  <  Rg  and  X2 

a.  \\\>*B 

The  solution  takes  the  form, 


la  + 


iT  »  exp  a^t,  where  a'X  +  (-0)2):K2and  a2  *  X  -  (-U)2)1^2 

which  is  an  egpenifitgiii.  function  whose  dominant  term  is  an  exponentially 
rising  term. 


l«*l 


<  R  and 
s 


>.2  >  ia  ♦  a 


& 


iT  =  A^  exp  a^t;  both  a^  and  a^  are  positive,  so  again  we  have  an 
exponentially  rising  current. 
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Summary  of  Analysis 


2  r  2 

Oscillations  occur  if  uu  is  positive  |^if  X  < 


C* 


1.  decreasing  oscillation  amplitude  if  X  <  0 

2.  Increasing  amplitude  if  X  >  0 

3.  constant  amplitude  if  X  »  0 

2 

Exponentially  varying  current  if  is  <  0 

1.  decreasing  exponential  if  X  <  0 

2.  increasing  exponential  if  X  >  0 

3.  decreasing  or  increasing,  depending  on  arbitrary  constants, 
thus  on  inital  conditions,  if  X  =  0 


where  X 


“ 1/2  6s 


37 


DISTRIBUTION 


Cn— finding  General 

U.  8.  Amy  Hate  riel  Command 

Washington,  D.  C. 

Attn:  AMCPP-CM-MI 

Office  of  the  Chief  of  Research  fc  Development 
Department  of  the  Army 
Washington  25,  D.  C. 

Commanding  Officer 
Pica tinny  Arsenal 
Dover,  New  Jersey 
Attn:  Tech  Library 

Commanding  Officer 

U.S.  Army  Research  Office  (Durham) 

Box  CM,  Duke  Station 

Durham,  North  Carolina  -  3  copies 

Commander 

U.S.  Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring  19,  Maryland 
Attn:  Tech  Library 

Department  of  the  Navy 
Bureau  of  Naval  Weapons 
Washington  25,  D.  C. 

Attn:  DLI-2,  Tech  Library 

Commander 

Naval  Research  Laboratory 
Washington,  D.  Co 
Attn:  Tech  Library 

DDC  Headquarters 
Defense  Documentation  Center 
Cameron  Station  Building  No.  5 
5010  Duke  Street 
Alexandria,  Virginia 

Attn:  TISIA  -  20  copies 

National  Bureau  of  Standards 
Washington,  D.  C. 

Attn:  Library 


39 


DISTRIBUTION  (Con't) 


University  of  Michigan 

Department  of  Electrical  Engineering 

Azin  Arbor,  Michigan 

Attn:  Professor  W.  G.  Dow 

University  of  California 
Department  of  Physics 
Berkeley,  California 

Attn:  Professor  L.  B.  Loeb 

Thomas  A.  Edison  Industries 
Edison  Laboratories 
West  Orange,  New  Jersey 
Attn:  R.  W.  Crowe 

Raytheon  Manufacturing  Company 
Receiving  Tube  Division 
55  Chapel  Street 
Newton  58,  Massachusetts 
Attn:  P.  W.  Stutsman 

University  of  Minnesota 
Dept  of  Physics 
Minneapolis,  Minnesota 

Attn:  Professor  H.  J.  Oskam 

Bell  Telephone  Laboratories,  Inc. 

Murray  Hill,  New  Jersey 

Attn:  A.  D.  White/M.  A.  Townsend 

Lockheed  Missile  and  Space  Co. 

3251  Hanover 
Palo  Alto,  California 
Attn:  Dr.  Leon  Fisher 

Minneapolis:  Honeywell  Regulator  Company 
500  Washington  Aye.,  South 
Hopkins,  Minnesota 
Attn:  B.  T.  McClure 

Sandia  Corporation 
Albuquerque,  New  Mexico 
Attn:  G.  W.  McClure 


DISTRIBUTION  (Don't) 


Microwave  Laboratory 
V.  W •  Hansen  Laboratories  of  Physics 
Stanford  University 
Stanford,  California 

Attn:  Dr„  F.  W.  Crawford 

General  Electric  Company:; 

Nela  Park  Laboratories 
Cleveland,  Ohio 

Attn:  Dr.  Carl  Kenty 

INTERNAL  DISTRIBUTION 

Horton,  B.  M./McEvoy,  R.  W.,  Lt  Col 

Apstein,  M./Gerwln,  H.  L„/Guarino,  P.  A./Kalmus,  H.  P. 

Schwenk,  C*  C. 

Hardin,  C.  D.,  Lab  100 
Sommer,  H.,  Lab  200 
Hatcher,  R«  D.,  Lab  300 
Hoff,  R.  S.,  Lab  400 
Nilson,  Ho,  500 

Flyer,  I.  N.,  Lab  600 

Campagna,  J.  H./Apolenls,  C.  J.,  Div  700 

DeMasi,  R. ,  Div,  800 

Landis,  P.  E.,  Lab  900 

Seaton,  J.  W.,  260 

Van  Trump,  J.  H.,  930 

Kaiser,  Q.  C.,  920 

Ward,  A,  L. ,  930  -  15  copies 

Reddan,  M.  J.,  930 

Belknap,  D.  J.,  930 

Hers  hall,  P.,  930 

Young,  R.  To,  320 

Harris,  F.  T.,  310 

Botkin,  I. /Godfrey,  T.  B./Eichberg,  R.  L. 

Bryant,  W.  T./Distad,  M  F./McCoskey,  R.  E  ./Moo  re  head,  J.G. 
Brenner,  M.,  723 

Technical  Information  Office,  010  -  5  copies 
HDL  Library  -  5  copies 

Qualified  requesters  may  obtain  copies  of  this  report  from  DDC 
5010  Duke  Street,  Alexandria,  Virginia 


